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1.  Foreword  (optional) 

The  aim  of  this  work  was  to  test  and  validate  a  new  approach  to  quantitatively  analyze  well 
hydrographs  in  shallow,  fractured  and  karstic  aquifers.  This  approach  was  demonstrated  to  be 
reliable,  and  estimates  of  transmissivity  (T)  and  specific  yields  (Sy)  in  different  portions  of  a  flow 
regime  can  be  obtained  by  the  use  of  this  less  costly,  non-invasive,  passive  method. 

Under  the  ARO  funded  project,  field  work  included  collecting  hydrograph  data  and 
conducting  aquifer  testing  at  selected  wells  at  three  sites:  the  Ammunition  Burning  Ground  in 
Crane,  IN;  the  Y-12  Plant  in  Oak  Ridge,  TN,  and  the  Main  Cantonment  area  of  Fort  Campbell, 
KY.  The  results  of  this  work  indicate  that  the  hydrograph  analysis  method  produces  aquifer 
parameter  estimates  (matrix  T)  of  comparable  reliability  to  those  obtained  using  conventional 
aquifer  testing  techniques.  Hydrograph  analysis  results  differ  from  those  of  slug  tests  in  cases 
where  cavities  are  located  in  the  well's  completion  zone.  This  results  because  the  slug  test  results 
are  biased  toward  the  higher  conductivity,  near  well-bore  conduits,  whereas  the  hydrograph 
method  evaluates  the  lower  conductivity,  matrix  T. 
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4.  Statement  of  the  Problem 

Numerous  Army  facilities  have  historically  conducted  activities  that  have  contaminated 
groundwater  with  both  organic  and  inorganic  materials.  In  order  to  properly  manage  these  sites, 
aquifer  parameters  must  be  obtained  for  adequate  characterization  and  predictive  modeling  to  be 
conducted.  Some  of  these  sites  are  located  in  karst  aquifers,  which  was  the  focus  of  this  work.  If 
the  technique  being  developed  here  could  be  demonstrated  to  be  reliable  at  these  karst  sites,  more 
reliable  aquifer  information  could  be  obtained  at  a  reduced  cost,  thus  decreasing  the  overall  cost 
of  characterizing  and  remediating  groundwaters  at  Army  facilities. 

Army  sites  are  not  the  only  areas  in  which  groundwater  has  been  contaminated  by  past  or 
current  activities.  Many  other  DoD  and  DOE  facilities,  as  well  as  privately  owned  properties, 
have  contaminated  groundwater  systems  requiring  characterization  and  remediation.  Many  of 
these  sites  are  also  located  in  karst  terrains,  and  hence,  the  hydrograph  analysis  method  could 
provide  more  reliable  data  at  a  reduced  cost  for  the  site  facilities  management  and  their 
consultants. 

Analytical  techniques  are  well  developed  for  hydrographs  obtained  from  streams  and  springs, 
where  data  are  cast  in  terms  of  total  discharge.  In  contrast,  well  hydrographs  are  plots  of  water 
level  versus  time.  It  was  hypothesized  that  three  segments  on  a  recession  curve  from  wells  in  a 
fractured  or  karst  aquifer  (multiple  porosity  systems)  represent  drainage  from  three  types  of 
storage:  conduit  or  larger  fractures,  smaller  fractures,  and  matrix.  A  method  was  developed  to 
estimate  matrix  transmissivity  (T)  and  specific  yields  (Sy)  in  these  different  components  of  the 
aquifer  using  well  hydrograph  data. 

Prior  to  this  funded  work,  the  well  hydrograph  analysis  method  had  only  been  applied  to  a 
limited  extent  at  one  site,  and  application  and  refinements  at  other  fractured  or  karst  sites  were 
needed  to  determine  if  the  method  could  be  more  universally  applicable  as  a  standard  analysis 
tool  in  other  shallow,  fractured  (multiple  porosity)  aquifers.  The  main  benefits  of  this  method 
are  in  the  acquisition  of  matrix  T  and  Sy  values  for  multiple  components  of  a  fractured  system 
rather  than  solely  for  the  continuum  as  is  determined  from  more  traditional  aquifer  tests.  In 
addition,  the  method  involves  passive  monitoring  of  a  well  bore,  which  is  desirable  in 
contaminated  areas  because  typical  well  testing  requires  pumping  water  out  of  the  borehole. 
Containment  and  treatment  of  the  pumped  waste  water  are  costly,  and  use  of  the  hydrograph 
analysis  method  results  in  avoidance  of  these  costs.  Hence,  if  the  hydrograph  method  could  be 
demonstrated  to  be  applicable  at  other  sites,  additional  information  on  aquifer  parameters  could 
be  obtained  at  a  reduced  cost. 

The  shape  of  the  rising  limb  of  a  hydrograph  is  largely  dictated  by  the  characteristics  of  a 
storm  event,  whereas  the  shape  of  the  recession  limb  is  largely  independent  of  the  character  of 
the  storm  (Linsley  et  al.,  1982).  Recession  limb  analysis  often  leads  to  two  or  more  line 
segments  that  represent  responses  in  the  different  portions  of  the  ground-water  system:  (1)  a  fast 
response  to  conduit  or  fracture  flow;  (2)  slower  responses  due  to  flow  through  smaller  fracture 
networks  and  unfractured  porous  media  (White,  1988).  At  relatively  shallow  depths  (<50  m)  in 
areas  with  sufficient  precipitation,  rapid  water  level  responses  are  expected  in  conduits  and 
fractures,  in  contrast  to  slower  responses  in  the  more  diffuse  portions  of  the  aquifer.  Because 
fractured  and  conduit  systems  have  little  hydraulic  resistance  in  contrast  to  porous  media, 
recharged  water  is  expected  to  drain  quickly.  The  porous  media  part  of  the  system  has  much 
lower  hydraulic  conductivity  and  responds  more  slowly  to  transient  events. 
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Relatively  few  studies  (Atkinson,  1977;  Rorabaugh,  1960)  have  yielded  quantitative  data  on 
aquifer  parameters  using  well  hydrographs.  It  has  been  previously  noted  that  stream  flow 
recession  curves  can  be  approximated  by  three  straight  lines  on  a  semilogarithmic  plot,  with  the 
lines  representing  three  different  types  of  storage:  stream  channels,  surface  soil,  and  groundwater 
(Barnes,  1940;  Linsley  et  al.,  1982).  Work  has  been  conducted  using  well  hydrographs  in  the 
karst,  Cambrian  Maynardville  Limestone  at  the  Y- 12  Plant  in  Oak  Ridge,  Tennessee.  Shevenell 
(1996)  hypothesized  that  three  segments  on  a  recession  curve  from  wells  in  a  karst  or  fractured 
aquifer  also  represent  three  types  of  storage:  larger  fractures  or  conduits,  smaller  fractures,  and 
matrix  portions  of  the  aquifer. 

The  purpose  of  the  ARO  funded  work  was  to  test,  validate  and  refine  the  scheme  to 
quantitatively  analyze  well  hydrographs  in  shallow,  fractured  rock  and  karstic  aquifers.  If  this 
scheme  could  be  demonstrated  to  be  reliable  at  other  sites,  estimates  of  transmissivity  and 
specific  yield  in  different  portions  of  the  flow  regime  could  be  obtained  by  use  of  an 
inexpensive,  non-invasive  method  that  does  not  disturb  the  aquifer. 


5.  Summary  of  the  Most  Important  Results 

Results  based  on  field  observations  and  theoretical  considerations  show  water  level 
variations  in  wells  in  response  to  storm  events  are  similar  to  those  monitored  in  springs 
discharging  in  karst  aquifers.  Methods  have  long  been  used  to  analyze  spring  hydrograph 
recession  curves,  and  the  work  presented  here  suggests  these  methods  can  be  extended  to  analyze 
the  similar  well  hydrograph  recessions.  In  this  work,  hydrograph  analysis  results,  which  provide 
an  estimation  of  T  in  matrix  intervals,  were  compared  to  results  obtained  from  traditional  aquifer 
testing  methods. 

Table  1  lists  all  of  the  wells  monitored  and  evaluated  as  part  of  the  funded  project.  The 
hydrograph  analysis  method  could  only  be  tested  against  a  relatively  small  number  of  results 
obtained  from  pumping  tests.  During  this  study,  no  wells  could  be  pumped  at  the  Crane  site, 
only  one  could  be  pumped  at  the  Ft.  Campbell  site,  yet  several  wells  at  the  Y- 12  site  could  be 
pumped.  However,  at  all  three  sites  numerous  wells  were  tested  with  the  use  of  hydrographs, 
which  do  not  require  withdrawal  of  ground  water,  and  this  is  precisely  one  of  the  main  reasons  to 
investigate  the  non-invasive  hydrograph  method. 

Table  2  summarizes  results  of  the  estimation  of  aquifer  T  at  all  three  sites.  Arithmetic 
means  were  computed  at  individual  wells  in  cases  where  there  were  multiple  T  estimates  for  the 
well.  Geometric  mean  values  were  computed  when  averaging  T  from  multiple  well  locations 
because  this  parameter  is  usually  considered  to  be  spatially  log-normally  distributed.  At  all  three 
sites,  results  of  the  hydrograph  analysis  technique  closely  match  (within  an  order  of  magnitude) 
those  of  traditional,  invasive  aquifer  parameter  estimation  techniques  (aquifer  pumping  tests  and 
slug  tests  conducted  in  matrix  intervals).  Slug  tests  conducted  in  wells  completed  in  cavities 
show  much  higher  estimated  T  values  as  a  result  of  the  larger  component  of  quick  flow  through 
conduits  relative  to  matrix  intervals  (Tables  2  and  3).  The  fact  that  differences  in  T  values  occur 
between  techniques  may  be  explained  by  the  difference  in  scale-of-measurement  between 
techniques.  For  example,  slug  testing  is  usually  considered  a  valid  indicator  test  for  T  within 
close  radial  proximity  to  the  tested  well,  whereas  pumping  tests  provide  a  T  estimate  for  a  larger 
portion  of  the  aquifer  measured  radially  outward  from  a  well.  Hydrograph  analysis  may  provide 
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aquifer  T  estimates  for  a  larger  area  than  that  of  pumping  or  slug  tests,  not  radially  from  a  well, 
but  upgradient  and  dependent  on  the  aquifer  drainage  geometry. 

The  results  of  this  work  support  the  hypothesis  that  well  hydrographs  may  be  used  to 
quantitatively  assess  the  hydraulic  properties  of  a  well-developed,  submerged,  fractured  or  karst 
aquifer.  Results  from  particular  wells  indicate  that  the  method  may  be  useful  in  areas  that  do  not 
contain  cavities,  but  that  do  have  multiple  porosities  that  are  drained  (e.g.,  matrix  plus  varying 
fracture  sizes/porosities).  Wells  in  larger  cavities  appear  to  influence  the  hydrograph  T  to  a 
higher  degree  than  purely  matrix  values  and  this  aspect  will  be  evaluated  as  part  of  future  work. 

There  are  limitations  to  this  methodology.  Sharp  storm  pulses  and,  hence,  well-defined 
hydrograph  recessions  with  multiple  limb  slopes  are  required  to  make  useful  quantifications. 
Complete  recessions  must  occur  before  the  hydrograph  is  influenced  by  the  next  storm. 
Simplifications  in  conduit  geometry  were  assumed  so  that  general  relationships  between  aquifer 
properties  could  be  obtained,  although  these  simplifying  assumptions  were  not  incorporated  into 
the  quantitative  analysis.  Instead,  quantitative  analysis  focuses  on  the  concept  of  karst  aquifer 
storage  depletion  as  a  whole  based  on  the  work  of  Rorabaugh  (1964),  and  conduit  drainage 
(independent  of  specific  conduit  shape)  as  recorded  by  the  well  hydrograph  during  periods  of 
storm  recession. 

The  aquifer  parameter  estimation  technique  to  obtain  matrix  T  presented  here  is  an 
alternative  to  the  commonly  used  pumping  and  slug  testing  methodologies.  This  method 
provides  realistic  estimates  of  aquifer  parameters  without  the  need  of  stressing  the  aquifer,  and 
hence  provides  information  during  times  of  natural  flow  conditions.  The  conduit,  fracture  and 
matrix  portions  of  the  aquifer  upgradient  of  the  monitoring  point  are  represented  in  this 
technique,  providing  a  better  understanding  of  aquifer  behavior  than  would  be  obtained  using 
point  measurement  techniques  alone.  A  similar  argument  can  be  made  for  the  case  of  discharge 
from  a  prominent  spring.  Such  hydrographs  provide  for  an  understanding  of  the  upgradient 
aquifer  a  spring  drains  by  separately  accounting  for  primary  and  secondary  water  storage  rather 
than  grouping  all  permeabilities  together.  The  work  presented  here  suggests  that  the  well 
hydrograph  analysis  technique  could  prove  useful  in  multiporosity  areas  where  pumping  test  data 
are  lacking  or  where  ground  waters  are  contaminated  and  pumping  is  either  too  costly  or  simply 
not  permitted. 

In  using  this  well  hydrograph  analysis  method,  commonly  collected  hydrograph  data  are 
used  to  go  beyond  conventional,  qualitative  methods,  and  begin  to  quantify  some  aspects  of  the 
karst  aquifer.  The  results  presented  here  from  49  wells  at  three  different  karst  sites  were 
compared  with  those  of  more  traditional  aquifer  testing  methods  (pumping,  slug).  What 
transmissivities  mean  in  any  highly  heterogeneous  aquifer  is  a  contentious.  However,  at  all  three 
sites,  the  hydrograph  T  (which  estimate  matrix  T)  agree  quite  well  with  data  obtained  from 
pumping  tests  from  the  same  wells,  or  from  slug  tests  in  nearby  matrix  dominated  wells  (Table 
3).  What  the  T  from  hydrograph,  pump  or  slug  test  results  represent  in  the  context  of  a  karst 
aquifer  is  debatable,  yet  the  latter  two  methods  are  frequently  used  now,  and  the  hydrograph  data 
provide  the  same  level  of  information  as  other,  more  traditional  methods  that  stress  the  aquifer. 
Use  of  the  well  hydrograph  technique  in  determining  matrix  T,  in  combination  with  (1) 
traditional  tracer  test  results  to  obtain  flow  velocities  in  quick  flow  portions  of  the  aquifer,  (2) 
slug  tests  in  quick  flow  dominated  portions  of  the  aquifer  to  estimate  T  of  the  conduit/fracture 
zones,  and  (3)  estimates  of  the  percentage  of  the  aquifer  to  which  to  apply  the  differing  Sy  values 
will  allow  improved  characterization  of  hydrologic  parameters  in  heterogeneous,  multiple 
porosity  systems  at  Army  and  many  other  karst  sites. 
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10.  Appendixes 


Table  1.  Summary  of  transmissivity  values  (m2/d)  computed  with  the  hydrograph  analysis  method,  slug  tests,  and  pumping  tests. 
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Table  2.  Comparison  of  transmissivities  estimated  by  different 
methods  for  individual  wells  at  the  three  sites. 


Well 

Total 

Zone 

Average 

Pumping 

Average 

Depth  (m) 

Type 

Hydrograph 

T  (m2/d) 

Slug 

T  (mVd) 

T  (m2/d) 

Crane  site 

03-32 

13.7 

C? 

0.1 

102 

03-34 

15.2 

C? 

1.5 

151 

03c24 

8.6 

C 

1.2 

11.6 

Average  of  4  nearby  matrix  wells  (Murphy,  1995): 

0.8 

Fort  Campbell  site 

2m5d 

28.6 

F 

0.09 

0.09 

0.1 

2mw4 

29.9 

F/C 

0.08 

0.3 

5mw3 

24.4 

C 

0.47 

31 

15mw5 

30.3 

c 

23 

455 

47mw3 

6.1 

c 

0.008 

581 

Oak  Ridge  site 

GW-056 

16.8 

9 

3.6 

2.1 

GW-057 

7.6 

? 

4.5 

1.3 

GW-603 

22.9 

M 

10.1 

7.6 

GW-684 

39.5 

c 

0.5 

0.6 

97.7 

GW-685 

42.2 

F/C 

3.2 

4.0 

GW-714 

44.2 

M 

8.7 

6.9 

GW-715 

13.6 

C 

21.5 

0.4 

GW-734 

31.4 

C 

22.4 

5.4 

50.8 

GW-735 

25.3 

C 

33.0 

0.3 

Average  of  5  nearby  matrix  wells: 

0.8 

Zone  Type:  M  -  matrix,  F  -  fracture,  C  -  cavity. 


Table  3.  Summary  of  computed  transmissivities  estimated  by  different  methods  at  the  three  sites. 


Average  Matrix  T  (m1 2/d) 

Average 

Slug 

Pumping 

Hydrograph 

Slug  Test 

Test 

Test 

Analysis 

T  (m2/d)  in  Cavities 

Crane  site 

8.4 

0.8* 

1.2 

32  ±63 

Number  of  wells  tested 

2 

4 

5 

6 

Total  number  of  tests 

16 

4 

45 

.  54 

Ft.  Campbell  site 

0.3 

0.1 

0.8 

51+226 

Number  of  wells  tested 

4 

1 

13 

9 

Total  number  of  tests 

24 

1 

48 

50 

Oak  Ridge  site 

l.l2 

2.2 

5.6 

70  ±33 

Number  of  wells  tested 

5 

8 

31 

2 

Total  number  of  tests 

30 

8 

64 

5 

Averages  are  geometric  mean  T  values  computed  for  all  wells  in  the  particular  category  (slug,  pumping, 
hydrograph).  The  total  number  of  wells  tested  for  each  type  of  test  (e.g.,  slug)  is  listed  on  the  first  line 
following  the  geometric  mean  T  value.  The  total  number  of  the  particular  test  (e.g.,  slug)  conducted  at  all 
wells  tested  is  listed  on  the  second  line  following  the  T  value. 

1  Summarized  from  Murphy,  1995. 

2  Summarized  from  Jones,  1997,  unpublished  data,  and  Shevenell  (unpublished  data). 


